Middle-Field Cusp Singularities in the Magnetization Process of One-Dimensional 

Quantum Antiferromagnets 
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We study the zero-temperature magnetization process {M—H curve) of one-dimensional quantum 
antiferromagnets using a variant of the density-matrix renormalization group method. For both the 
5 = 1/2 zig-zag spin ladder and the S = 1 bilinear-biquadratic chain, we find clear cusp-type 
singularities in the middle-field region of the M — H curve. These singularities are successfully 
explained in terms of the double- minimum shape of the energy dispersion of the low-lying excitations. 
For the 5=1/2 zig-zag spin ladder, we find that the cusp formation accompanies the Formi- liquid 
to non-Fermi-liquid transition. 

PACS numbers: 75.10.Jm, 75.40.Cx, 75.30.Cr 



Low-dimensional antiferromagnetic (AF) quantum 
spin systems with various spin magnitude S and various 
spatial structures, have been a field of active researches, 
both experimentally and theoretically. In particular, the 
magnetization process (M — H curve, M: magnetization, 
H: magnetic field) of AF spin chain at low-temperatures 
have recently drawn much attention, because it exhibits 
various phase-transition-like behaviors: e.g., the criti- 
cal phenomena AM ~ ^/H — associated with the 
gapped excitation (excitation gap oc He) or with 

the saturated magnetization (s.t the saturation field H^), 
magnetization plateau, |^ and, the first-order tran- 
sition. ^ These are field-induced phase transitions of the 
ground state, which reflect non-trivial energy-level struc- 
ture at zero field. 

What we consider in this Letter is another type of 
singularity which has not been discussed so much: the 
cusp singularity dX H — H^usp in the middle-field re- 
gion {He < i?cusp < ^^s)- Existence this type of sin- 
gularity was first demonstrated by Parkinson for 
the integrable Uimin-Lai-Sutherland (or, SU(3)) chain, 
and has also been known for some other integrable 
spin chains, [p"2|-p^ and ladders. Derivation of the 

magnetization cusp for each model, however, relies on 
the model's integrability in an essential manner, restrict- 
ing the Hamiltonian to be of somewhat unrealistic form. 
Hence, whether or not this type of singularity can be 
found in realistic systems is a highly non-trivial ques- 
tion. In the present Letter, we make the first system- 
atic numerical study of the "middle-field cusp singular- 
ity" (MFCS, for short) in the zero-temperature M — H 
curve for non-integrable systems, through which we give 
a positive answer for the above question. Namely, we 
show that the S — 1/2 zig-zag spin ladder and the S' = 1 
bilinear-biquadratic chain actually exhibit MFCS in the 
M — H curve. 

The numerical method we employ is the product- 
wavefunction renormalization group (PWFRG), p^ ] 
which is a variant of the density-matrix renormalization 



group (DMRG). ^ Efiiciency of the PWFRG in cal- 
culation of the M ^ H curve has been demonstrated in 
previous studies; [|j20) the PWFRG wiU allow us to ob- 
tain the M — H curve in the thermodynamic limit with 
enough accuracy to detect "weak" (non-divergent) singu- 
larity hke the MFCS. 

Consider the 5 = 1/2 zig-zag spin ladder, whose ac- 
tual realization can be made as a quasi-one-dimensional 
material. The M — H curve of the system with bond 
alternation has recently been studied, where the cusp sin- 
gularity has not been discussed. |Q The Hamiltonian of 
the system is given by 
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[Si ■ Si+i + JSi ■ Si+2] H^'^ S^ , (1) 



where Si is the S = 1/2 spin operator at the i-th site. We 
have normalized the nearest neighbor coupling to unity, 
and have denoted the next-nearest coupling by J (> 0). 
In Fig. 1, we show the M — H curve calculated by us- 
ing the PWFRG. We see the M - H curve for J > 1/4 
exhibits a clear MFCS. 

We take down-spin-particle picture, where one down 
spin in the saturated (all up) state is regarded as a par- 
ticle. Then the system near Hg can be regarded as that 
of interacting Bose particles, which reduces to the well- 
known delta-function Bose gas (5-BG) model |23| in the 
low-energy limit. Important point is that the S- 

BG at low density is equivalent to spinless free-Fermi 
gas. [Ullljl^ Hence, for discussion of the M — H curve 
near Hg, we may treat the system as that of the spin- 
less Fermions which is completely characterized by the 
one-particle excitation energy dispersion uj{k). 

The one-down-spin excitation energy a>(fc) is calculated 
to be 



(jj(fc) = cos fc — 1 -t- J(cos 2k — 1), 



(2) 



which we depicted in Fig. 2. It should be noted that, at 
J = 1/4, there occurs a qualitative change in the shape of 



1 



Lu{k): For J < 1/4, u}{k) have a single minimum at fc = tt, 
while, for J > 1/4, uj{k) at fc = tt changes into a local 
maximum and two minima newly appear (correspond- 
ing fc-positions are determined from cos A; = —1/(4 J)). 
| p6| , p4[ Then, the van Hove singularity associated with 
the double-minimum shape of uj{k) gives a simple expla- 
nation of the MFCS in the M - H curve for J > 1/4. 


Let us make a quantitative analysis as follows. From 
the bottom of uj{k), we obtain Hg — 2 for 1/4 > J > 0, 
and Hs = I + 2J + 1/(8 J) for J > 1/4. As the apphed 
field H is decreased below Hs , the down-spin density be- 
comes to be non-zero. The M ~ H curve is obtained 
from 



H = 



1/2- — / R{k)dk, 
27r J 


(3) 


1- J cj{k)R{k)dk, 


(4) 


dE{M) 
dM ' 


(5) 



where R{k) is the zero-temperature Fermi distribution 
function which is unity inside the "Fermi vacuum" but is 
zero otherwise. Hence, how the particles fills the energy 
band tells us the essential behavior of the M — H curve. 

In < J < 1/4, the M — H curve is smooth in the 
whole range ofO < M < 1/2. The one-particle energy 
Lu{k) near the bottom position k — tt has the expansion 
uj{k) = -2+ i(l - 4J)(AA:)2 + ^(16J - l)(Afc)4 + ... 
where we have introduced Ak = k — ir. The H — M curve 
near is then calculated from (||), (Q), and (H) to be 



Hs-H 



7r2(l -4J) 



(AM)2 -I- 



7r4(16J- 1) 
24 



{AM f ■ 



(6) 



where AM = 1/2 — M . In the correctly mapped 5-BG 
treatment for a class of AF spin chains, ^ there may be 
(AAf)^ term due to the finiteness of the effective cou- 
pling. Hence, for the present model, we fit the PWFRG- 
calculated M — H curve with 

Hs-H = a{AMf [1 + 7(AAf ) + 6{AMf . . .] , (7) 

to check whether the obtained value of a agrees with the 
free-Fermion prediction. The best-fit results are a = 2.6 
(J ~ 0.1) and a = 0.92 (J = 0.2), which are consistent 
with the free-Fermion prediction: a = 2.960 . . . ( J = 0.1) 
and a = 0.9869 .. . (J = 0.2). 

At J = 1/4, (AA:)^-term in uj{k) vanishes, lead- 
ing to a different form of the expansion: Hs — H = 

The best-fit value 



a{AMy 



1 + ^{AM) + S{AM)^ ■ 
of a from the PWFRG result lies in the range a = 
14 ~ 17, 1 25 1 which is in reasonable agreement with the 
free-Fermion prediction a — 12.2.... Hence the free- 
Fermion picture also holds for J = 1/4, and we have 



AM ^ {H — HsY^^ with the critical exponent 1/4 being 
different from the "standard" value 1/2, [ p|~^ J^ support- 
ing the finite-size diagonalization result. [ ^6[ 

We have thus seen that the system for J < 1/4 is 
a Fermi liquid, near Hg. Hence, also for J > 1/4, we 
may expect that the system continues to behave as a 
Fermi liquid. Qualitatively, the Fermi-liquid character 
well explains the M — H curve, in particular, the appear- 
ance of the MFCS. For quantitative discussion, however, 
there emerges an important difference from the case of 
J < 1/4: Due to the double- minimum shape of a;(fc), the 
system becomes to be two-component liquid [each com- 
ponent is composed of modes around each minimum]. 
Since the separation into two components may not be 
complete, there should remain interactions between the 
components. Such a correlated multicomponent system 
may often behave as a non-Fermi liquid, or, Tomonaga- 
Luttinger (TL, for short) liquid, [p7|-p9| whose typical 
example is the Hubbard chain. |^ 

The TL liquid is characterized by the smooth edge 
of the momentum distribution at the Fermi point kp- 
R{k) - (fcp - fc)'^ (fc ~ /fcp, C > ). Hi Then, for very 
small particle density, we assume 



R{k) = Roikp - \k\)'^/k- 



(8) 



where Rq is a constant, and we have assumed that i?(0) 
remains to be finite even in the vanishing particle den- 
sity {kp — > 0). Assuming that each component has a 
parabolic energy dispersion uj{k) = aAk^ {Ak = k — k*) 
around each minimum k = k* , we have the following form 
of the ground-state energy density Eq as a function of 
the total particle number density p: 



Eg 
Ctl 



Ctl^P , 
Ctl(-Ro, C) 



24(C -1-1)2 



i?^(C + 2)(C + 3) 



(9) 
(10) 



For the two-component (non-interacting) Fermi liquid, 
we have C = 1 and Rq — 2 giving Ctl(-Ro7 0) = 1. There- 
fore, deviation of Ctl from unity implies the non-Fermi 
liquid character of the system. 

Since p corresponds to AM = 1/2 — A/, the TL- 
liquid expression (^ gives the coefficient a in (^ as 
a = CTLCT7r2/4 where a is calculated from the band cur- 
vature at each minimum of tjj{k). Explicitly, we have 
cr = (16 J2 - 1)/(8J) and 



(16 J2 - 1) 2 
32J ■ 



(11) 



The best-fit values of a from the PWFRG calculation 
are 5.9 ± 0.3 (J = 0.4) and 7.4 ± 0.4 (J = 0.5). J2|] 
These values disagree with the free-Fermion values ((pT[) 
with Ctl = 1): a = 1.20 . . . ( J = 0.4) and a = 1.85 . . . 
(J = 0.5). Rather, our calculation leads to Ctl ~ 4 in 



2 



(^T|), showing a clear sign of the non-Fermi-hquid char- 
acter of the system. Remarkably, in the Hubbard chain 
with on-site interaction U (> 0) which can be viewed as 
a two-chain S = 1/2 spin ladder via the Jordan- Wigner 
transformation, the same factor Ctl — 4 appears in the 
U — oo limit ||3l[] or in the low-density limit. ||3^ There- 
fore, we may conclude that, near H^, the S = 1/2 zig-zag 
ladder for J > 1/4 is a two-component non- Fermi TL liq- 
uid. 

Let us next consider the bilinear-biquadratic (BLBQ, 
for short) chain with the Hamiltonian 



'^BLBQ — ^2 \_^^ 



Si+1 + [i{Si 



hY^SI (12) 



where Si — {Sf, Sf , Sf) is the S = 1 spin operator at the 
site i. 11 In Fig. 3 we show PWFRG-calculated M-H 
curves for (3 = 0.45, 0.6, 0.8, 1.0, whare clear MFCS is 
seen. Q 

In Ref. 1^] , we made a quantitative test for the square- 
root critical behavior, ||l|,|,|| M - A{H - H^)^/'^ [A: 
amplitude) , of the M — H curve near the lower critical 
field He (which is proportional to the excitation gap at 
H = 0). There, we found a critical value (3c (~ 0.41) 
where the critical exponent changes from 1/2 to 1/4, 
which is in close similarity to the J = 1/4 case of the 

= 1/2 zig-zag spin ladder. We should remark that the 
change in the shape of w(fc) at (3c has recently been found. 
Hi Therefore, the MFCS in the BLBQ chain is weU ex- 
plained in terms of the Fcrmi/Tomonaga-Luttinger- liquid 
picture. We should also note that the observed two- 
component Fermi/Tomonaga-Luttinger- liquid behavior 
for (3 > (3c, 10 is also consistent with the appearance 
of the MFCS. 

To summarize, in this Letter we have studied the 
middle-field cusp singularity (MFCS) in the zero- 
temperature magnetization process {M — H curve) for 
antiferromagnetic spin systems in one dimension. For 
the S* = 1/2 zig-zag spin ladder and the S — 1 bilinear- 
biquadratic chain, we have found clear MFCS in the 
M—H curve obtained by using the product- wavefunction 
renormalization group method which is a variant of the 
density-matrix renormalization group (DMRG) method. 
We have explained the mechanism for the MFCS in terms 
of the shape-change in the energy dispersion curve of the 
low-lying excitation. Further, for the 5 = 1/2 zig-zag 
spin ladder, we have shown that the formation of the 
MFCS accompanies the Fermi-liquid to non-Fermi-liquid 
(Tomonaga-Luttinger-liquid) transition in the character 
of the system. 

As far as we know, what we have found for the S = 
1/2 zig-zag spin ladder is the first non-trivial example 
of physically observable MFCS. For actual experi- 
mental observation, we should of course take the finite- 
temperature effect into account, which can be made by 
the "finite-temperature DMRG" . p9(| As an important 



implication drawn from the present study, we should note 
that the essential mechanism for appearance of the MFCS 
is the multi-minimum structure of the low-lying excita- 
tion energy. Such structure may be well expected for 
systems with non-trivial spatial structures and/or com- 
peting interactions, which often accompany incommen- 
surability in physical quantities. In fact, for the S = 1/2 
zig-zag spin ladder in the large J region (J > 0.5) (in- 
commensurability at zero field is reported ^^), we have 
observed another cusp near the lower critical field, whose 
details will be published elsewhere. 
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FIG. 1. The M — H curve of the zig-zag spin chain cal- 
culated by the PWFRG with number of the retained bases 
m = 30. (a) J = 0.1, 0.2, and 0.25. (b) J = 0.35, 0.4, and 
0.5. Inset: magnification of the curves around the cusps. 
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FIG. 2. The dispersion curve of the "one-down-spin" for 
the zig-zag spin chain. 
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FIG. 3. The M - H curve of the BLBQ chain for (3 = 0.45, 
0.6, 0.8 and 1.0 calculated by the PWFRG with number of 
the retained bases m — 70. 
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